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Stablization  of  Nanotwinned  Microstructures  in  Stainless  Steels  Through  Alloying  and 

Microstructural  Design 

Raymundo  Arroyave,  Texas  A&M  University 

1  Preamble 

One  of  the  most  significant  challenges  for  the  development  of  high-performance  structural 
materials  is  the  de-stabilization  of  optimized  microstructures  when  operating  at  relatively  high 
homologous  temperatures.  While  this  problem  is  especially  pressing  in  areas  of  power  generation 
and  aerospace  propulsion,  the  problem  of  stabilization  of  nano-structures  is  highly  relevant  to 
many  other  applications  of  interest  to  ARO.  In  this  STIR  we  focused  on  the  generation — through 
thermo-mechanical  processing —  and  stabilization — through  alloying  and  nano-precipitate 
engineering — of  nanotwinned  stainless  steels.  This  preliminary  work  has  provided  valuable 
insight  into  the  mechanisms  responsible  for  the  stabilization  of  a  wider  class  of  nanotwinned 
materials.  Knowledge  gained  from  this  program  is  currently  being  used  to  develop  guidelines  for 
the  design  of  stable  nanotwinned  microstructures. 

The  microstructure  design  will  exploit  nano-twinning-induced  grain  boundary 
engineeringstrategies  which  will  target  high  density  of  low  energy  grain  and  nano-twin 
boundaries,  pinned  by  nano-scale  precipitates  -  formed  at  very  high  temperatures-  for  high- 
temperature  strength  and  stability.  The  uniqueness  of  the  proposed  approach  originates  from  the 
synergistic  combination  of  computational  alloy  design  expertise  with  microstructural  design  and 
thermo-mechanical  processing  capabilities  for  accelerated  discovery  of  the  next  generation 
advanced  high  temperature  stainless  steels. 

2  Research  Activities 

In  this  STIR  we  proposed  to  investigate  the  microstructure  and  alloy  composition  on  the 
formation  and  stability  of  nanotwinned  microstructures.  The  long-term  alloy  design  strategy 
focuses  on  microalloying  additions  to  1)  reduce  the  stacking  fault  energy  (SFE),  enhance 
twinning  formation,  and  increase  twin  volume  fraction;  2)  reduce  the  twin  size  via  increasing 
interstitial  content  while  also  taking  into  account  solid  solution  hardening  in  twin  size;  3)  form 
high  temperature,  mainly  intergranular,  carbides,  nitrides,  carbonitrides  to  pin  twin  boundaries, 
increasing  their  high  temperature  stability;  4)  form  nano-scale  intermetallic  precipitates  in  the 
presence  of  Co  and  refractory  elements  for  creep  resistance  and  high  temperature  strength. 
Microstructural  design  strategy  would  concentrate  on  twinning-induced  grain  boundary 
engineering  and  distribution  of  nano-precipitates/carbides.  The  work  would  also  involve  a 
fundamental  study  on  the  thermal  stability  of  deformation  nano-twins  at  high  temperatures  and 
their  influence  on  the  recovery/recrystallization  of  austenitic  stainless  steels. 

In  the  present  STIR  project  focused  mainly  on  the  formation  and  stability  of  deformation 
nano-twins  in  the  single  crystals  of  low  SFE  austenitic  stainless  steels  ,  experimentally.  We  did 
some  preliminary  investigation  on  the  phase  stability  of  precipitate  phases  likely  to  stabilize 
twinned  structure  at  elevated  temperatures.  This  computational  work  was  accompanied  by  few 
experimental  validation  cases  through  microstructural  investigations  and  thermo-mechanical 
testing.  The  experimental  work  was  performed  in  collaboration  with  Prof.  Ibrahim  Karaman. 
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2.1  Fundamental  Investigations  on  Twinning  Behavior  in  Fe-based  Alloy  Systems 

The  likely  operative  deformation 


FCC 


TWIN 


o 


B 


Figure  1:  Schematic  diagram  of  twining  in  fee 
systems. 


mechanism  of  the  steel  alloys  can  be 
understood  in  part  by  examining  the  behavior 
stacking  fault  energy  (SFE),  which  is  generally 
defined  as  the  deformation  characteristics  and 
mechanical  behavior  of  the  materials.  In  a 
simple-minded  way,  from  low  to  high  SFE,  the 
plastic  deformation  of  steel  alloys  can  be 
divided  into  three  preferred  stages:  (1) 
martensitic  transformation,  (2)  twinning  and 
gliding  and  (3)  gliding  [1-3].  This  means,  that 
with  lower  SFE,  the  alloy  possesses  higher 
tendency  to  twin  in  FCC  phase.  Although  there 
have  been  a  great  number  of  experimental  and  theoretical  investigations  on  the  intrinsic  factors 
that  control  the  SFE,  many  questions  remain.  This  is  particularly  true  when  considering  that  one 
of  the  major  goals  of  the  present  project  were  to  incorporate  a  metric  for  twinnability  in  the 
overall  alloy  design  scheme. 

From  a  thermodynamic  point  of  view,  twinning  in  fee  systems  corresponds  to  a  local  phase 
transition  where  the  normal  ABC  stacking  is  locally  transformed  to  a  AB  stacking  corresponding 
to  an  hep  crystal  structure.  Thus,  to  a  first  approximation,  SFE  can  in  principle  be  expressed  as  a 
function  of  molar  surface  density,  molar  Gibbs  free  energy  change  during  the  phase 
transformation  (in  this  case,  fee  to  hep),  and  surface  energy  of  the  interface  (  y  =  2pAG  +  2a). 
Thermodynamics  models  have  been  built  to  estimate  the  required  free  energy  change  for  the 
local  fee  to  hep  transformation  (essentially  a  martensitic  transition).  From  these  models,  it  has 
been  concluded  that  by  adding  Ni,  A1  and  Cu,  SFE  goes  higher;  whereas  by  adding  Cr,  SFE 
decreases  [3,  5].  This  indicates  the  importance  of  the  chemical  composition. 

Another  important  issue  is  the  criterion  to  increase  the  twinning  density.  Twinnability  is  the 
quantity  characterizing  the  ability  of  a  material  to  deform  through  twinning  deformation  [6].  This 
quantity  can  also  be  expressed  as  a  function  of  the  unstable  stacking  energy  (USE)  and  unstable 
twinning  energy  (UTE).  From  first  principle  calculations,  these  two  properties  can  be  evaluated. 

This  model  provides  the  qualitatively  results 
comparing  to  the  experiments. 

Thermodynamic  Approach 
The  stacking  fault  energy  (SFE)  in  austenitic 
stainless  steel  can  be  estimated  using  either 
thermodynamic  models  or  first-principle 
calculations.  In  this  work  the  thermodynamic 
database,  TCFE6  V6.2,  is  utilized  to  calculated  the 
SFE  with  the  model  [4,  7] 

T  =  2.88xlO“5AG^£  +  16  (mj  ■  moF1) 
where  16  mJ-moF 1  represents  the  interfacial 
wjw.%  °  °  energy  of  the  twinning.  To  evaluate  this  model,  we 

Figure  2:  The  SFE  calculated  by  the  attempted  to  compare  the  experimental  results  on 
thermodynamic  model  comparing  to  the  Fe22Mno.6C  based  alloy  with  various  Cr  content 
experimental  results[4] 
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performed  by  Dumay  and  collaborators  [4]  with  the  results  obtained  from  the  thermodynamic 
approach.  As  shown  in  Fig.  2  there  is  a  discrepancy  between  the  calculation  using  TCFE6  (V6.2) 
and  the  published  results.  Our  model  calculations  overestimate  the  SFE  in  the  alloy  system  by  a 
factor  of  two.  While  the  reason  for  the  discrepancy  is  not  clear  at  the  moment,  it  is  likely  that  the 
major  contributor  is  the  fact  that  the  free  energy  differences  between  fee  and  hep  according  to  the 
thermodynamic  database  are  overestimated.  This  is  likely  to  the  fact  that  the  Fe-rich  hep  phase 
has  not  been  fully  assessed.  This  is  in  part  due  to  the  scarce  thermodynamic  and  phase  equilibria 
data  where  hcp-Fe  is  a  relevant  phase.  Given  the  fact  that  this  quantity  (difference  in  energy 
between  fee  and  hep)  is  not  easily  accessible  experimentally,  we  proceeded  to  investigate  the  use 
of  DFT  calculations  to  estimate  SFE. 


Density  Functional  Theory  Calculations 

The  fcc-hcp  martensitic  transformation  is  presented  from  atomic  point  of  view  in  Fig.  3.  The 
initial  atomic  arrangement  is  in  sequence  of  ABCABCA  and  transformed  into  ABCACAB  with 
respect  to  the  inhabit  plane  of  atom  A  [8].  It  also  represents  the  model  for  calculating  the 
intrinsic  stacking  fault  as  a  function  of  the  difference  in  energy  between  the  faulted  and  perfect 
crystal  structure,  and  the  specific  surface  area  of  the  plan  parallel  to  the  fault  [9]: 

Fsf  ~  F0 
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‘2D 
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(.)  -  •  •  • 


[m]  I 


Figure  3:  The  schematic  process  of  the  martensitic  transformation  [8] 


While  the  calculation  of  the  quantity  above  can  be  routinely  done  using  first  principles 
approaches,  other  approximations  can  be  made  in  order  to  facilitate  the  rapid  estimation  of  SFE 
within  minimum  computational  effort.  According  to  the  Axial  Next-Nearest-Neighbour  Ising 
(ANNNI)  model,  the  total  energy  of  a  system  with  different  sacking  sequences  can  be 
represented  as  [9]: 


E  =  E0 


SA 


i+3 


The  energy  different  of  the  stacking  fault  can  therefore  be  estimated  based  on  the  free 
energies  of  the  HCP  and  FCC  phases  in  a  first  order  approximation  [10,  11] 

Fsf  ~  Fo  =  4/i  —  0(J2 )  ~  Fhcp  —  Ffcc 

This  model  has  been  utilized  to  study  the  effects  of  substitutional  [10,  12,  13]  and  interstitial 
[8,  14]  elements  to  SFE  of  the  austenitic  stainless  steels.  Also,  the  finite  temperature  model  is 
implemented  for  practical  applications  [15,  16].  While  this  method  is  rather  trivial  to  implement 
in  conventional  DFT  codes,  the  fact  that  alloys  can  only  be  simulated  through  the  use  of 
supercells  makes  the  calculation  of  the  difference  in  energy  as  a  function  of  alloy  composition 
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problematic.  In  order  to  overcome  this  limitation,  we  have  examined  two  alloys,  Fe74.5Cri3.5Nii2 
and  Fe7o.5Cri7.5Nii2,  using  the  EMTO-CPA  method[17].  Fig.  4  shows  the  optimum  wigner-seitz 
radii  of  FCC  and  HCP  phases  of  these  two  alloys  and  the  SFEs  are  calculated  accordingly  as 
shown  in  Fig.  4.  Comparing  to  the  work  by  Vitos  et  al.  [11],  the  calculations  present  reasonable 
results. 


Figure  4:  The  stability  of  the  FCC  and  HCP  phases  of  Fe74.5Cri3.5Nii2  and  Fe70.5Cr17.5M12 


Table  1.  The  calculated  energy  states  and  SFE;  n  stands  the  magnetic  moment  (average  of  Fe  and 
Mn)  of  FCC  and  HCP  (HbY  Yo  -  Ym.rt  -  Y  are  the  stacking  fault  energies  (mJ/m2)  at  0  Kelvin,  affected  by 
magnetic  moment  contributions  to  entropy  at  room  temperature  and  total  SFE  at  room  temperature. 


YLfcc 

^HCP 

Yo 

Ym,rt 

Y 

Fe74.5crl3.5M12 

1.61 

0.00 

7.03 

30.71 

37.75 

[14] 

1.62 

0.00 

8.50 

36.20 

44.60 

Fe70.5Cr17.5M12 

1.54 

0.00 

6.67 

29.84 

36.51 

[14] 

1.54 

0.00 

-1.10 

29.70 

28.60 

Fe65.5Cri7.5Nii2Mn5 

1.32 

0.00 

7.19 

26.97 

34.15 

The  results  shown  in  Table  1  are  worth  further  examination:  as  can  be  seen  from  the  table, 
the  OK  stacking  fault  energy  calculated  at  OK  is  rather  law.  Moreover,  it  can  be  seen  how 
different  the  magnetization  is  between  fee  and  hep.  This  means  that  magnetic  properties  are 
likely  to  play  a  major  role  due  to  their  contribution  to  the  entropy  differences — and  free  energy 
differences.  In  fact,  the  table  shows  that  magnetic  contributions  can  account  for  more  than  75% 
of  the  stacking  fault  energy  in  Fe-based  compositions  close  to  those  of  typical  stainless  steels. 

In  order  to  further  investigate  the  twinning  deformation  behavior  in  Fe-based  fee  alloys,  we 
proceeded  to  investigate  the  Generalized  Stacking  Fault  Energy  surface  according  to  the  method 
proposed  by  Janhatek  and  collaborators[18].  The  methodology  allows  for  the  investigation  of 
shearing  deformation  parallel  to  the  (1 1 1)  planes  along  different  important  directions,  [11-2]  and 
[-110]  (Fig.  5a).  Two  different  types  of  deformations — affine  and  alias — have  been  investigated. 
Under  an  affine  shear  deformation,  all  atoms  are  shifted  parallel  to  the  shearing  direction  by  a 
distance  proportional  to  their  distance  from  the  fixed  basal  plane.  In  the  alias  regime,  only  the  top 
layer  is  displaced  in  the  shearing  direction.  Fig.  5b  contains  a  schematic  of  the  deformation 
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mechanism  considered  in  this  work. 


Figure  5a:  Geometry  of  shear  deformation  on 
(111)  planes[18]. 


Figure  5b:  The  schematic  view  of  possible 
deformation  regimes:  (a)  affine  shear,  (b)  the  alias 
shear[18]. 


Figure  6:  Generalized  stacking  fault  energy  surface  for  pure  Fe,  without  considering 
magnetism.  Green  curve  corresponds  to  fcc-hcp  transformation  under  alias  shear. 
Calculations  were  done  at  OK. 

Fig  .  6  shows  the  Generalized  Stacking  Fault  Energy  (GFSE)  surface  for  Fe  along  the  [11-2] 
and  [-110]  directions  under  alias  and  affine  shears.  The  curve  corresponding  to  the  alias  shear 
along  the  [11-2]  direction  corresponds  to  the  conventional  deformation  driving  the 
transformation  from  fee  to  hep.  This  curve  shows  that  the  SFE  is  actually  negative  at  the  absolute 
minimum  of  the  transformation  path.  This  point  corresponds  to  the  Intrinsic  Stacking  Fault.  The 
negative  value  means  that  the  structure  is  unstable  against  this  deformation.  This  makes  sense  as 
the  hep  structure  in  Fe  when  no  spin  magnetization  is  considered  is  actually  more  stable  than  fee 
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Figure  7:  Generalized  stacking  fault  energy  (GSFE)  surface  for  pure  Fe  along  the  [11-2]  direction 
considering  different  magnetic  states  at  OK. 

Fig.  7  shows  similar  calculations  as  shown  in  Fig.  6,  comparing  different  magnetic  states  in 
fee  Fe:  NM  (non  magnetic),  FM  (anti-ferromagnetic)  and  AFM  (anti-ferromagnetic).  Here  it  is 
important  to  note  that  the  actual  magnetic  state  for  fee  Fe  is  more  complex  than  the 
configurations  considered.  However,  it  is  important  to  note  that  regardless  of  the  state,  hep  is  still 
more  stable  than  fee.  Note,  however,  that  in  the  case  of  AFM  state,  the  barrier  necessary  to 
overcome  in  order  for  the  transformation  to  occur  is  actually  higher  than  in  the  NM  state.  Note 
also  how  a  FM  arrangement  leads  to  a  less  stable  hep.  Note  also  that  in  this  case  no  contributions 
due  to  entropy  have  been  taken  into  account.  It  is  likely  that  these  effects  are  sufficient  to  raise 
the  GFSE  curve,  although  it  is  unlikely  that  the  ISF  reaches  positive  values  under  the 
approximations  considered. 

Future  work  along  the  lines  outlined  above  will  consist  of  including  all  possible 
contributions  to  the  free  energy  of  the  GFSE  surface.  Moreover,  alloying  effects  will  be 
considered. 

2.2  Plastic  Deformation  Mechanisms  in  316  Stainless  Steels 

To  help  the  alloy  design  by  providing  better  understanding  of  microstructural  evolutions 
through  microstructural  characterization,  single  crystals  of  stainless  steel  316  were  fabricated  to 
study  the  effect  of  mechanical  deformation  on  twin  formation  and  the  changes  in  the  twin 
structure  during  recovery  and  recrystallization.  Single  crystals  were  grown  using  the  Bridgman 
technique  in  He  and  were  then  cut  into  dog-bone  shape  tensile  samples  with  the  gage  dimensions 
of  8  mm  x  3mm  x  1mm  using  wire  electro  discharge  machining. 

316  single  crystals  with  composition  shown  in  Table  2  were  used  for  tension  test  at  room 
temperature  and  -80°C.  The  elemental  composition  table  was  provided  using  EDS  analysis  and, 
therefore,  may  include  a  ±5  wt%  error.  Nitrogen  was  added  to  this  composition  to  study  its  effect 
on  twin-ability  of  stainless  steel  316L  and  forming  possible  nitro-carbides. 


“NM 

“FM 

“AFM 


Relative  Displacement  along  [11-2] 
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Table  2.  Composition  of  316LN  single  crystals  from  EDS  analysis 


Element 

C 

N 

Cr 

Cu 

Mn 

Mo 

Ni 

O 

P 

S 

Si 

Fe 

wt% 

0.08 

0.20 

17.73 

0.25 

1.07 

2.26 

11.78 

0.0012 

0.014 

0.014 

0.44 

bal. 

316  [111]  single  crystals  were  cut  into  dog  bone-shaped  specimens  and  were  heat  treated  for  12 
hours  at  1080°C,  then  at  1100°C  for  1  hour  and  consequently  water  quenched,  in  order  to 
homogenize  the  microstructure.  One  sample  was  then  used  in  uniaxial  tension  test  for  20% 
deformation  at  room  temperature  (Fig.  8a).  EBSD  results  of  this  sample  are  shown  in  Fig.  8b. 
The  results  show  that  a  twinning  texture  of  {111}/{1 12}  exists,  and  the  band  of  the  twins  vary 
from  less  than  1pm  to  larger  than  40pm. 


Figure  8:  (a)  Strain-Stress  behavior  of  3 16  [111]  crystal  sample  uniaxial  tension  test  at  room  temperature 
(b)  EBSD  analysis  of  3 16  [1 1 1]  single  crystals  deformed  to  20%  at  room  temperature. 


Figure  9:  TEM  images  of  20%  deformed  3 16  [1 1 1]  single  crystal  samples 
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The  20%  deformed  316  [111]  single  crystal  samples  was  studied  using  TEM.  Fig.  9  shows 
two  TEM  images  of  this  sample.  In  Fig.  9-a,  a  lOnm  wide  band  can  be  detected,  and  Fig.  9-b 
shows  the  diffraction,  indicating  that  this  is  twinning.  Also  Error!  Reference  source  not 
found.9-c  shows  another  set  of  bands,  around  lOOnm  wide,  with  the  corresponding  diffraction 
pattern  in  Fig.  9-d  it  is  suggested  that  these  bands  may  be  stacking  faults.  These  TEM  and  EBSD 
results  are  promising  and  investigation  will  continue  on  optimizing  the  twin  size  and  its 
interaction  with  precipitates,  as  well  as  other  mechanical  properties  at  room  temperature  and  high 
temperature 

2.3  Phase  Stability  in  Stainless  Steels 

As  mentioned  above,  the  formation  of  nano-size  precipitates  is  beneficial  to  stabilize 
nanotwinned  microstructure  at  elevated  temperatures.  It  is  recently  shown  in  316  stainless  and 
Fe-Mn  austenitic  steels  that  nano-carbides  prefer  forming  at  grain  boundaries  and  faults  [19-21] 
and  that  these  phase  can  stabilize  nanostructures  and  improve  properties  at  elevated 
temperatures [22,  23].  However,  one  problem  is  to  control  the  temperature  at  which  these 
carbides  form.  Low  temperature  carbides/nitrides  should  be  suppressed  such  that  the 
microstructure  is  stable  during  long  high  temperature  exposure.  To  achieve  this,  we  have 
considered  using  elements  such  as  Si  (to  a  certain  extent  not  to  form  sigma  phase)  and  Al,  and/or 
forming  strong  substitutional  -  interstitial  couples  with  elements  such  as  Mn.  Formation  of  nano¬ 
precipitates  of  MX  (M=Nb  or  W,  X  =  N  or  C)  at  very  high  temperatures  can  suppress  the 
precipitation  of  other  phases  at  lower  temperatures  [24],  enhancing  the  high  temperature  stability 

and  creep  resistance  of  austenitic 
stainless  steels  at  desired  operating 
temperatures.  Addition  of  W,  Ta, 
and  Nb  should  also  provide 
significant  solid-solution 

strengthening  in  the  proposed 
steels. 


1000  1500 

T,  Kelvin 
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Figure  10:  Phase  analysis  for  Stainless  Steel  316  with  Based  on  the  composition  for 

composition  shown  in  Table  1.  Stainless  Steel  316  (shown  in 

Table  1),  the  equilibrium  phases  are  analyzed  and  shown  in  Fig.  10,  according  to  the  TCFE6 
V6.2  database,  where  Sigma  phase  is  the  FeCr  inter-metallic  compound.  It  indicates  that  the 
Laves  phase  is  not  stable  in  this  alloy  at  temperatures  higher  than  600K.  In  order  to  suppress  the 
formation  of  M23C6,  Nb  and  Ti  are  added  to  316  ,  as  alloy  2  and  3,  to  form  Fe(Nb,  Ti)C.  The 
effect  of  temperature  and  weight  percent  of  Ti  and  W  on  the  formation  of  ferrite  (BCC), 
austenite  (FCC),  Laves  phases,  and  NbTi  carbide  is  shown  in  Fig.  11. 


Alloy  2 


Alloy  3 


Ti  case  V  case  W  case 


(b) 
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Figure  11:  The  effect  of  Ti  and  W  additions  on  the  volume  fraction  of  (a) 
Austenite  (FCC),  (b)  Laves  Phase,  and  (c)  TiNbC. 


It  can  be  observed  that  the  Ti  is  most  effective  to  stabilize  Austenite  and  TiNbC  but  it  is  not 
as  important  as  V  to  Laves  phase.  These  contour  plots  can  be  utilized  for  selecting  the  alloy  and, 
more  importantly,  for  defining  the  domain  of  the  material  optimization.  In  future  work,  we  will 
incorporate  other  elements  and  we  will  examine  the  stability  of  other  possible  stabilizing  nano¬ 
phases,  including  intermetallic  compounds  (such  as  aluminides). 

3  Outlook 

The  ultimate  goal  of  this  preliminary  project  was  to  understand  the  effects  of  alloying  and 
microstructure  design — through  thermo-mechanical  processing — on  the  size,  density  and 
stability  of  nanotwins  in  stainless  steels.  This  STIR  forms  the  basis  of  a  longer-term  program 
aimed  at  the  design  of  new  austenitic  stainless  steels  with  ultrahigh  strength,  ductility  and — as  a 
result  of  the  stable  nanotwinned  microstructure — high-temperature  strength.  The  long-term  goal 
of  this  effort  is  to  develop  austenitic  stainless  steels  capable  of  operating  at  temperatures  up  to 
850°C  that  demonstrate  yield  strength  levels  on  the  order  of  more  than  a  GPa,  ultimate  tensile 
strengths  (UTS)  of  more  than  2  GPa,  and  specific  strength  of  more  250  kNm/kg  at  ambient 
temperature. 
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